Supersonic Air Flow due to Solid-Liquid Impact 
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A solid object impacting on liquid creates a liquid jet due to the collapse of the impact cavity. 
Using visualization experiments with smoke particles and multiscale simulations we show that in 
addition a high-speed air-jet is pushed out of the cavity. Despite an impact velocity of only 1 m/s, 
this air-jet attains supersonic speeds already when the cavity is slightly larger than 1 mm in diameter. 
The structure of the air flow resembles closely that of compressible flow through a nozzle - with the 
key difference that here the "nozzle" is a liquid cavity shrinking rapidly in time. 



PACS numbers: 47.55.D-, 47.60.Kz, 47.11. St, 47.80.Jk 
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Taking a stone and throwing it onto the quiescent sur- 
face of a lake triggers a spectacular series of events which 
has been the subject of scientists' interest for more than 
a century [ll-[l7j : upon impact a thin sheet of liquid (the 
"crown splash") is thrown upwards along the rim of the 
impacting object while below the water surface a large 
cavity forms in the wake of the impactor. Due to the 
hydrostatic pressure of the surrounding liquid this cavity 
immediately starts to collapse and eventually closes in a 
single point ejecting a thin, almost needle-like liquid jet. 
Just prior to the ejection of the liquid jet the cavity pos- 
sesses a characteristic elongated "hourglass" shape with 
a large radius at its bottom, a thin neck region in the 
center, and a widening exit towards the atmosphere. 

This shape is very reminiscent of the converging- 
diverging ("de Laval") nozzles known from aerodynam- 
ics as the paradigmatic picture of compressible gas flow 
through, e.g., supersonic jet engines. In this Letter we 
use a combination of experiments and numerical simula- 
tions to show that in addition to the very similar shape, 
also the structure of the air flow through the impact cav- 
ity resembles closely the high-speed flow of gas through 
such a nozzle. Not only is the flow to a good approx- 
imation one-dimensional, but it even attains supersonic 
velocities. Nevertheless, the pressure inside the cavity is 
merely 2% higher than the surrounding atmosphere. The 
key difference, however, is that in our case the "nozzle" is 
a liquid cavity whose shape is evolving rapidly in time - 
a situation for which no equivalent exists in the scientific 
or engineering literature. 

Our experimental setup consists of a thin circular disc 
with radius Rq = 2 cm which is pulled through the 
liquid surface by a linear motor mounted at the bot- 
tom of a large water tank [l6| with a constant speed 
of Vb = 1 m/s. To visualize the air flow we use small 
glycerin droplets (diameter roughly 3 /im) produced by 
a commercially available smoke machine (skytec) com- 
monly used for light effects in theaters and discotheques. 
Before the start of the experiment the atmosphere above 
the water surface is filled with this smoke which is conse- 
quently entrained into the cavity by the impacting disc. 



A laser sheet (Larisis Magnum II, 1500mW) shining in 
from above illuminates a vertical plane containing the 
axis of symmetry of the system. A high-speed camera 
(Photron SAl.l) records the motion of the smoke parti- 
cles at up to 15,000 frames per second. Cross-correlation 
of subsequent images allows us to extract the velocity of 
the smoke which faithfully reflects the actual air speed 
[l8|. Our setup obeys axisymmetry and we use cylindri- 
cal coordinates with z = the level of the undisturbed 
free surface. 

In the beginning of the process (see the snapshot in 
Fig.[T](a)) air is drawn into the expanding cavity behind 
the impacting object with velocities of the order of the 
impact speed. In a later stage however, this downward 
flux is overcompensated by the overall shrinking of the 
cavity volume resulting in a net flux out of the cavity. 
The cavity shape at the moment when the flow through 
the neck reverses its direction is illustrated in Fig. [T](b). 
Towards the end of the cavity collapse a thin and fast 
air stream is pushed out through the cavity neck which 
is illustrated in Fig. [1] (c). From images such as those 
in Fig. [U we can directly measure the air speed u up to 
about 10 m/s as is shown in the inset of Fig. O 

In order to determine the flow speed at even higher 
velocities we revert to multiscale numerical simulations. 
Our numerical method proceeds in two stages: an incom- 
pressible stage at the beginning and a compressible stage 
towards the end of the impact process. During the first 
stage both air and liquid are treated as incompressible, 
irrotational, and inviscid potential fluids. To solve for the 
flow field and to calculate the motion of the interface we 
use a boundary integral method (BIM) as described in 
jl6| with extensions to include the gas phase [22]. At the 
moment that the air flow through the neck reverses, see 
Fig. [1] (b), the simulation enters into the second, com- 
pressible stage: from now on only the liquid motion is 
computed by the incompressible BIM. 

To simulate the air flow in the second stage we need 
to take compressibility into account meaning that a sim- 
ple potential flow description is no longer possible. For- 
tunately, at the end of the incompressible stage the air 
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FIG. 1: (a) After the impact of the disc an axisymmetric 
cavity is formed in its wake and air is entrained into this 
cavity, (b) Due to hydrostatic pressure from the surrounding 
Uquid the cavity starts to cohapse and the air flow reverses 
its direction, (c) As the collapse proceeds air is pushed out 
of the shrinking cavity at very high speeds. In (a)-(c) we 
overlaid images of the cavity shape (recorded with backlight) 
and images of the smoke particles (recorded with the laser 
sheet and artificially colored in orange). In the latter, the 
area illuminated by the vertical laser sheet is restricted by 
the minimum cavity radius ^T§\. A corresponding movie can 
be found in [Tol ]. 
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FIG. 2: The speed of the gas flowing through the neck (red 
curve) as a function of the shrinking cavity neck taken from 
the fully compressible simulations. The main plot demon- 
strates that sonic speeds are attained with the cavity pres- 
sure (blue curve) being less than 2% higher than the atmo- 
spheric pressure. The enlargement (inset) shows that the 
numerical scheme (red curve) agrees very well with the ex- 
perimentally measured velocity (black diamonds; the hole in 
the data between rneck = 16 mm and 22 mm is due to mea- 
surement uncertainties at low absolute velocities [3). Slight 
non-axisymmetric perturbations 20, 21] in the experimental 
setup may be responsible for the somewhat slower air speed 
of the experiment as compared to the simulation. One can 
clearly see the inversion of the flow direction from negative 
(into the cavity) to positive (out of the cavity) velocities. 



velocity profile is almost perfectly one-dimensional along 
the axis of symmetry. We can therefore describe the gas 
dynamics by the ID compressible Euler equations [23| 
in analogy to gas flowing through a converging-diverging 
nozzle. In the Euler equations we include two additional 
terms accounting for the variation of the nozzle radius 
in time and space h^. For the numerical solution we 
use a Roe scheme [23, which is highly appreciated 
for its computational efficiency and ability to accurately 
capture shock fronts. 

The two-way coupling between the gas and the liquid 
domains is accomplished via (i) the interfacial shape and 
its instantaneous velocity which is provided by the BIM 
and serves as an input into the gas solver and (ii) the 
pressure which is obtained from the solution of the Eu- 
ler equations and serves as a boundary condition for the 
BIM. Above the location of the initial free surface the 
surface pressure of the BIM remains atmospheric. 

Combining our experiments with these numerical sim- 
ulations leads to the main result of this Letter contained 
in Fig. [2j the collapsing liquid cavity acts as a rapidly de- 
forming nozzle, so violent that the air which is pushed out 
through the neck attains supersonic velocities (red line). 
Our simulations show that the pressure inside the cavity 
which is driving this flow is less than 1.02 atmospheres 
(blue line). From the inset one can tell that our simu- 
lations are in good agreement with the smoke measure- 
ments over the entire experimentally accessible range. It 
is interesting to note that even towards the end of the 
process (when sonic velocities are reached) there is a net 
flux of air upwards through the cavity. If the process was 
governed merely by the collapse of the neck itself one 
would expect the air to be pushed out of the neck region 
in both vertical directions. This net flow thus underlines 
the important role of the dynamics of the entire cavity. 

To determine more precisely at what point the air flow 
through the neck becomes sonic we show in Fig. [3] (a) the 
evolution of the local Mach number, Maneck = '^neck/c 
(with the gas velocity u^eck and the speed of sound c), 
for discs impacting at 1 and 2 m/s. We find that the 
speed of sound is attained at cavity radii as large as 0.5 
mm for the lower impact velocity and 1.2 mm for the 
higher impact velocity. 

In a steady state one could expect from the (com- 
pressible) Bernoulli equation that these very high air 
speeds would cause a greatly diminished air pressure in 
the neck region. Despite the unsteadiness of our sit- 
uation, the data presented in Fig. [3] (b) indeed shows 
that the pressure Pneck decreases significantly once the 
neck has shrunk to a diameter of roughly 4 mm (for the 
1 m/s impact) while before that point it is practically 
atmospheric throughout. Classical steady-state theory 
[25[ for a converging-diverging nozzle predicts that when 
Maneck = 1 the pressure at the neck reaches a minimum 
value of 
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FIG. 3: (a) The evolution of the local Mach number at the 
cavity neck for different impact speeds (red: 1 m/s, blue: 2 
m/s). For the 2 m/s impact speed sonic flow is attained at 
a cavity radius of 1.2 mm. (b) The pressure at the neck di- 
minishes due to Bernoulli suction as the neck radius shrinks 
and air is forced to flow faster and faster. The minimum 
pressure lies at about 0.6pa which is attained when the Mach 
number reaches unity, (c) The experimental image shows a 
pronounced kink at the neck which is not captured by the 
smoothly rounded curve predicted by the simulation without 
air (cyan line) . Only the inclusion of air effects into the simu- 
lations (red line) is able to reproduce the kinked shape caused 
by the low air pressure at the neck as well as the shape of the 
cavity above the neck. 



with Pa the atmospheric pressure and 7 = 1.4 the isen- 
tropic exponent. As shown in Fig. [3] (b) our situation - 
although highly unsteady - exhibits a similar behavior 
with Pneck ~ 0.6pa cis the Mach number becomes of order 
unity. 

In Fig.[3](c) we illustrate how this low pressure gives us 
a handle to observe the consequences of the supersonic air 
speed in our experiments: despite the air being three or- 
ders of magnitude less dense than water, it is able to exert 
a significant influence even on the shape of the liquid cav- 
ity provided that its speed is high enough [26|, [l^ • From 
the experimental image it is clear that the free surface 
close to collapse no longer possesses a smoothly rounded 
shape but instead shows a significant increase in curva- 
ture at the minimum (a "kink"). While this feature is 
not present in a simulation neglecting the influence of air 
as those in [l6|, the inclusion of air effects allows us to 
capture quite accurately the cavity shape observed ex- 
perimentally. This gives strong evidence that in the ex- 
periment the air indeed becomes as fast as predicted by 
the simulations and produces a Bernoulli suction effect 
strong enough to deform the cavity. 

The positive sign of iXneck (see Fig. [2]) indicates that 
the gas flow is directed upwards at the neck. At the 
same time, the air at the bottom of the cavity is pulled 
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FIG. 4: (a) The vertical air velocity normalized by the local 
speed of sound Ma = u/c as a function of the vertical position 
(the corresponding cavity image is shown in the middle) for 
^neck = 0.9 mm: the profile exhibits a sharp peak approxi- 
mately at the height of the neck, (b) A close-up of the zone 
around the neck illustrates the steepening of the velocity pro- 
files towards pinch-off (numbers 1-5 correspond to neck radii 
between 0.9 mm (number 1, bright red) and 0.5 mm (num- 
ber 5, dark brown)) and the development of the shock front 
at roughly 0.1 ms before pinch-off. The neck position 2;neck 
corresponding to curve 5 is shown by the dashed line, (c) A 
close-up of the area below the neck shows the location of the 
gas flow stagnation point ^^stag (dashed line). 



downwards by the moving disc. An interesting conse- 
quence of this competition between cavity expansion at 
the bottom and cavity shrinking in the neck is the exis- 
tence of a stagnation point with = as can readily be 
observed in Fig. |4] (a) and its magnification in Fig. |4] (c). 

As can be seen in the inset of Fig. [5l the distance 
between the neck and the stagnation point is no larger 
than roughly 5 mm prior to cavity closure. Nevertheless, 
the pressure at the stagnation point equals the overall 
pressure inside the cavity which is roughly atmospheric 
during the whole process (see Fig. [2]). Recalling that 
Pneck ^ 0.6pa this rcsults in a tremendous vertical pres- 
sure gradient which of course affects the dynamics of the 
cavity wall: the flow of air is so strong that it can drag the 
liquid along resulting in an upward motion of the cavity 
neck just before the final collapse. That this effect is in- 
deed present in the simulations can be seen from the red 
line in Fig. [5l For comparison, the cyan curve demon- 
strates that a single fluid simulation neglecting the air 
dynamics would predict a monotonously decreasing po- 
sition. The experimental data however is in quantitative 
agreement with the compressible simulations. Together 
with the cavity shape shown in Fig. [3] (c) these results 
constitute an impressive - albeit indirect - demonstration 
of the credibility of our numerical predictions despite the 
fact that, understandably, it is not possible to directly 
measure (super-)sonic air speeds with our smoke setup. 
Furthermore they show that the perfectly axisymmetric 
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FIG. 5: The vertical position of the cavity neck relative to the 
final closure height Zc as a function of the shrinking neck ra- 
dius from experiment (black diamonds) , simulations with (red 
line) and without (cyan line) air dynamics. The experimen- 
tal data is in quantitative agreement with the compressible 
simulations, while clearly the simulation neglecting air fails 
to capture the upward motion of the minimum induced by 
the large pressure gradient between the stagnation point and 
the cavity neck. Experimental error bars are determined by 
the number of vertically neighboring pixels all sharing the 
same minimum radius. The inset shows the approach of the 
stagnation point to the neck. 



approach of the simulations is justified and, therefore, 
that supersonic gas velocities are reached before instabil- 



ities [20I, HH inevitably destroy the axisymmetry of the 
system. 

Looking more closely at the velocity profile above the 
neck (see Fig. |l](b)) one finds that it possesses a discon- 
tinuous jump: the signature of a shock front developing 
in the air stream. While such a shock front is a com- 
mon phenomenon in steady supersonic flows, here we are 
able to illustrate its development even in our highly un- 
steady situation when the gas velocity passes from sub- 
to supersonic. 

In conclusion, we showed that the air flow inside the 
impact cavity formed by a solid object hitting a liquid 
surface attains supersonic velocities. We found that the 
very high air speeds can be reached even though the pres- 
sure inside the cavity is merely 2% higher than the sur- 
rounding atmosphere. This is due to the highly unsteady 
gas flow created by the rapidly deforming cavity. We il- 
lustrated how the air affects the cavity shape close to 
the final collapse in two different ways: (i) the initially 
smoothly curved neck shape acquires a kink which can 
be attributed to a Bernoulli suction effect and (ii) the 
initially downward motion of the neck reverses its direc- 
tion and starts to travel upwards. The quantitatively 
consistent observation of both effects in numerics and 
experiment makes us confident that our rather involved 
numerical procedure truthfully reflects reality. 
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